Dynamic Critical Phenomena in Trapped Bose Gases 



V.I. Yukalov 1 , E.P. Yukalova 2 and V.S. Bagnato 3 

1 Bogolubov Laboratory of Theoretical Physics 
Joint Institute for Nuclear Research, Dubna 141980, Russia 

. 2 Department of Computational Physics 

, Laboratory of Informational Technologies 

kjPj! Joint Institute for Nuclear Research, Dubna 141980, Russia 

3 Instituto de Fisica de Sao Carlos, Universidade de Sao Paulo 



< 

o 



oo 
o 



Caixa Postal 369, Sao Carlos, Sao Paulo 13560-970, Brazil 



ABSTRACT 



> 
OO 

O . 

I/"") . Nonlinear dynamics of a trapped Bose-Einstein condensate, subject to the action of a resonant external 
field, is studied. This field produces a spatio-temporal modulation of the trapping potential with the 
frequency close to the transition frequency between the ground state and a higher energy level. The 
evolution equations of fractional populations display a kind of critical phenomena at a critical line on 
the manifold of the system parameters. It is demonstrated that there exists a direct analogy between 
dynamical instability at this line and critical phenomena at a critical line of a related averaged system. 
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A dilute cloud of trapped Bose-condensed atoms, at low temperatures, can be described by a wave 
function satisfying the Gross-Pitaevskii equation [1,2]. This equation is nonlinear due to the atomic 
interactions through a local delta potential. The mathematical structure of the equation is that of the 
nonlinear Schrodinger equation. Stationary states of the latter, because of the confinement caused by 
a trapping potential, are restricted to discrete energy levels. These stationary eigenstates form a set of 
solutions that can be called nonlinear modes, in analogy to linear modes that are solutions of a linear 
Schrodinger equation. The nonlinear modes can also be called coherent modes, since the wave function of 
the Gross-Pitaevskii equation corresponds to a coherent state of Bose-condensed atoms. One can term 
as well these modes as topological modes, emphasizing that the wave functions related to different energy 
levels have different spatial topology. 

The modal structure of the confined condensate states is in close analogy with nonlinear optical 
modes of nonlinear waveguide equations in optics [3,4]. And the methods of creating nonlinear coherent 
modes of trapped Bose atoms [5,6] are also similar to those employed in optics, where one uses specially 
prepared initial conditions or invokes an action of external fields. Dynamics of fractional populations, 
characterizing the occupation of coherent modes of Bose-condensed trapped atoms, displays as well 
many phenomena equivalent to those known in optics, for instance collapses, revivals, and Rabi-type 
oscillations. The nonlinear dynamics of processes coupling bosonic coherent modes have been studied in 
several papers considering the general case [5,6], antisymmetric modes [7,8], dipole topological modes in 
a two-component condensate [9,10], and vortex modes [11-14]. 
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The aim of the present paper is twofold: First, we demonstrate the existence of a critical line for 
the population dynamics of Bose condensates, which is due to the nonlinearity of evolution equations. 
Second, we show that there is an intimate relation between the studied nonlinear dynamical system and 
a related averaged system, so that the instabilities happening at a critical line of the dynamical system 
are the counterparts of critical phenomena occurring in the averaged system. 

The time-dependent Gross-Pitaevskii equation has the form 

<p , H(<p) = - A- V 2 + U(f ) + ,4iVM 2 , (1) 

where U(r) is a trapping potential; A = Airfi 2 a s /mQ- } a s is a scattering length; mo is mass; and N is 
the number of particles; V is a potential of external fields. The wave function ip is normalized to unity, 
\\(p\ \ = 1. The nonlinear topological modes are defined as the eigenfunctions of the nonlinear Hamiltonian, 
that is, they are solutions to the eigenproblem 

H(ip n )(p n (r) = E n ip n (r). 

It is worth stressing here the principal difference between the topological coherent modes and el- 
ementary collective excitations. The former are self-consistent atomic states defined by the nonlinear 
Gross-Pitaevskii equation. While the latter are the elementary excitations corresponding to small devi- 
ations from a given atomic state and described by the linear Bogolubov-De Gennes equations [1,2]. 

We assume that at the initial time the system was condensed to the ground-state level with an energy 
Eq. So that the initial condition to eq. (1) is ip(r, 0) = tpo(r). Suppose, we wish to couple the ground 
state (fo with another state ipj having a higher energy Ej. The best way for doing this, is clearly, by 
switching on an external field, say V = V(f) cos cot, oscillating with a frequency u) which is close to the 
transition frequency 

Uj = (Ej - E )/h, 

so that the detuning Ao; = to — uj be small, \Auj/uj\ <C 1. Then one can look for a solution to eq. (1) in 
the form 

<p(r, t) = c (t)Mr) e~ iEot / h + Cj {t) Vj {r) e~ iE ^ h . (2) 

The considered situation is analogous to the description of nonlinear resonant processes in optics [3,4]. 
The validity of the quasiresonant two-level approximation (2) for the time-dependent Gross-Pitaevskii 
equation has been confirmed from the general point of view [5,6] and also proved for several similar cases 
by direct numerical simulations of the Gross-Pitaevskii equation [8,9,11-13], the agreement between the 
two-level picture and the simulations being excellent. 

The coefficients Cj(t) define the fractional level populations rij(t) = |cj(t)| 2 . The equations for these 
coefficients can be obtained by substituting the presentation (2) in eq. (1). It is again worth noticing 
the similarity of this presentation with the slowly-varying-amplitude approximation commonly used in 
optics, if one treats the factors Cj(i) as slow functions of time, as compared to the exponentials in eq. (2), 
which implies that \dci/dt\ <C -Ej. This approximation, being complimented by the averaging technique 
[15], permits one to slightly simplify the evolution equations for Ci(t). In this way [5,6], we come to the 
equations 

— = -ia rijCo - -f3cj e 1 ^ 1 , -± = -tat n cj - - (3 c e , (3) 
in which the transition amplitudes 

ay ee A j J |^(F)| 2 (2|^(r)| 2 - |^(r-)| 2 ) dr , P=\J Vd(Wf)w(f) 



H(<p) + V 
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caused by the nonlinearity and by the modulating field, respectively, are introduced, and the abbreviated 
notation a = aoj, with setting, for simplicity, aoj = ctjo, is used. Note that the transition amplitude (3 
is nonzero only if the potential V(r) depends on space variables. The initial conditions to eqs. (3) are 
c (0) = 1 and c,(0) = 0. 

Usually, one solves such evolution equations with fixed parameters a, (3 and Aw, keeping in mind a 
particular experimental realization. Instead of this, we have studied the behavior of solutions to eqs. (3) 
in a wide range of varying parameters. It turned out that this behaviour is surprisingly rich exhibiting 
new interesting effects. 

First of all, it is easy to notice that the number of free parameters in eqs. (3) can be reduced to 
two by the appropriate scaling. For this purpose, we measure time in units of a" 1 and introduce the 
dimensionless parameters 

b=\P\/a, 5 = Auj/a. 
It is also evident that eqs. (3) are invariant with respect to the inversion 

a — ► —a, (3 — ► — (3, Auj — > — Aw, t — > —t. 

Therefore it is possible to fix the sign of one of the parameters, say a > 0, since the opposite case 
can be obtained by the inversion. For concreteness, we shall also keep in mind that (3 is positive. 
The dimensionless detuning is assumed to always be small, \8\ <C 1. And the dimensionless transition 
amplitude b is varied in the region < b < 1. We have accomplished a careful analysis by numerically 
solving eqs. (3). When parameter b is small, the fractional populations oscillate reminding the Rabi 
oscillations in optics, where \(3\ would play the role of the Rabi frequency. The amplitude of these 
oscillations increases with increasing b. It would be more correct to say that in our case there exist 
nonlinear Rabi oscillations, as far as eqs. (3) differ from the corresponding equations for optical two-level 
systems by the presence of the nonlinearity due to interatomic interactions. This nonlinearity not only 
slightly modifies the Rabi-type oscillations of the fractional populations but, for a particular relation 
between parameters, can lead to dramatic effects. By accurately analyzing the behaviour of solutions to 
eqs. (3), with gradually varying parameters, we have found out that there exists the critical line 

6 + 5-0.5, 

at which the system dynamics experiences sharp changes [6]. For example, if the parameter b = 0.4999 
is kept fixed, and the critical line is crossed by varying the detuning 8, the following phenomena occur. 
When the detuning is zero, the fractional populations display the Rabi-type oscillations. After slightly 
shifting the detuning to 8 = 0.0001 the top of nj{t) and the bottom of no(t) become flat, and the oscillation 
period is approximately doubled. A tiny further variation of the detuning to 8 = 0.0001001 yields again 
drastic changes, so that the upward cusps of nj(t) and the downward cusps of no(t) arise. The following 
small increase of the detuning to 8 = 0.00011 squeezes the oscillation period twice. Then, making 8 
larger does not result in essential qualitative changes of the population behaviour. All dramatic changes 
in dynamics occur in a tiny vicinity of the critical line. The same phenomena happen when crossing the 
line b + 8 ~ 0.5 at other values of parameters or if 8 is fixed but b is varied. 

The unusual behaviour of the fractional populations is due to the nonlinearity of the evolution equa- 
tions (3). Systems of nonlinear differential equations, as is known, can possess qualitatively different 
solutions for parameters differing by infinitesimally small values. The transfer from one type of solutions 
to another type, in the theory of dynamical systems, is, generally, termed bifurcation. At a bifurcation 
line, dynamical system is structurally unstable. 
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The second aim of our paper is to show that the found instability in the considered dynamical system 
is analogous to a phase transition in a statistical system. To elucidate this analogy for the present case, 
we have to consider the time-averaged features of the dynamical system given by eqs. (3). To this end, 
we need, first, to define an effective Hamiltonian generating the evolution equations (3). This can be 
done by transforming these equations to the Hamiltonian form 

idcQ jolt = dH e j f I 8cq , ialcj /alt = dH e f j / dc* , 

with the effective Hamiltonian 

H eff = a norij + ±(j3 e iA ^ Cj + 0* e'^c*^) . (4) 

An effective energy of the system can be defined as a time average of the effective Hamiltonian (4). For 
this purpose, eqs. (3) can be treated by means of the averaging technique [15], as it is described in detail 
in Ref. [5], which provides the guiding-center solutions. Substituting the latter in eq. (4), together with 
the time-averaged fractional populations, results in the effective energy 

E eff = (ab 2 /2e 2 )(b 2 /2e 2 + 5), 

where e is a dimensionless average frequency defined by the equation 

£ 4 (e 2_ 62) = (e 2_ 6 2_ e2(5) 2_ 

This effective energy, represents the time average of the effective Hamiltonian (4) and, thus, characterizes 
the average features of the system. As an order parameter for this averaged system, one can take the 
difference of the time-averaged populations, 

rj = no — n j = 1 — b 2 /e 2 . 

The capacity of the system to store the energy pumped in by the resonant field can be described by the 
pumping capacity 

Cp = dE eff /d\(3\. 

The influence of the detuning on the order parameter is characterized by the detuning susceptibility 

Xs = \drj/dd\. 

Analyzing the behaviour of the introduced characteristics as functions of the parameters b and 5, we 
found out that they exhibit critical phenomena at the critical line b + 5 = 0.5, which coincides with 
the bifurcation line for the dynamical system. Expanding these characteristics over the small relative 
deviation r = \b — b c \/b c from the critical point b c = 0.5 — 6, we obtain 

n- Vc c^(l-26y/ 2 , Cpc^r-V 2 , Xs c±r-V 2 , (5) 

where n c = rf(b c ) and r — > 0. As is seen, the pumping capacity and detuning susceptibility display 
divergence at the critical point. The related critical indices for Cp, n, and \& are equal to 1/2. These 
indices satisfy the known scaling relation: 

ind(Cp) + 2 ind(rj) + ind(xs) = 2, 



4 



where ind is the evident abbreviation for index. 

In order to clarify what is the origin of the found critical effects for the studied dynamical system, let 
us return back to eqs. (3). Again we pass to dimensionless notation measuring time in units of a. By 
means of the substitution 

C ° = (~~2^) 6XP {* ( q ° + \ *) } ' Cj = (~T^~) 6XP ( qi ~ 5 *) } ' ^ 
where p, qo, and q\ are real functions of t, equations (3) can be reduced to the form 

dp , H 2 • dq bp 

- = -b^l-p 2 S ,nq, -= P+7 ==co S q + 5, (7) 

in which q = q\ — q$. Note that this reduction to an autonomous dynamical system is valid for arbitrary 
detuning 5. Moreover, this system possesses the integral of motion 

I(p, q) = -p 2 — b\J 1 — p 2 cos q + Sp, 

which can be defined by using the initial conditions p(0) = — 1 and q(0) = corresponding to the 
conditions co(0) = 1 and Cj(0) = 0. This gives 

/(-l,0) = l-<5. 

The existence of the integral of motion means that the dynamical system is integrable in quadratures. 
This fact does not help much for studying the time evolution of the system, since the formal solutions p(t) 
and q(t) are expressed through rather complicated integrals, so that the system evolution, anyway, is to 
be analysed numerically. However, the property of integrability implies that the appearance of chaos in 
the system is impossible. Consequently, the observed critical effects in no way could be related to chaos. 
Then what is their origin? The answer to this question comes from the analysis of the phase portrait for 
eqs. (7) in the rectangle defined by the inequalities 

-l<p<l, 0<q<27r. 

This analysis shows that, if b + S < 0.5, then the motion starting at the initial point p(0) = —1, q(0) = 
is oscillatory, with a trajectory lying always in the lower part of the phase rectangle, below the separatrix 
given by the equation 

-p 2 — b\J 1 — p 2 cos q + bp — (3 = 0. 

When b + 5 = 0.5, the separatrix touches the initial point, so that the following motion occurs in the 
phase region above the separatrix. In this way, if we consider, under the given initial conditions, the 
parametric manifold formed by the parameters b £ [0, 1] and 5 < 1, then the critical line b + 5 = 0.5 
separates the parametric regions related to two different types of solutions to eqs. (7). 

Summarizing, we have considered the population dynamics of a trapped Bose-Einstein condensate, 
subject to the action of a resonant spatio-temporal modulation of the trapping potential. The consid- 
eration has been reduced to a two-level picture corresponding to the quasiresonant approximation. A 
careful analysis of the derived two-level evolution equations has been made for the wide range of vary- 
ing parameters. Such a variation of parameters can be easily realized by changing trap characteristics, 
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varying the number and kind of condensed atoms, and by changing the scattering lengths using Feshbach 
resonances [1,2]. It turned out that on the manifold of possible parameters there exists a critical line, 
where the evolution equations display structural instability. We have demonstrated that this instability 
for a dynamical system is analogous to a phase transition for a stationary averaged system. For the 
latter, one can define an order parameter, pumping capacity, and detuning susceptibility which exhibit 
critical phenomena at the critical line. The origin of this critical line is elucidated by showing that it 
divides the parametric manifold onto two regions corresponding to qualitatively different solutions of the 
evolution equations. 

The consideration of the nonlinear dynamics of trapped Bose-condensed gases can be based on the 
two-level approximation because of the resonant character of the modulating external field. However, 
the accuracy of this approximation remains undetermined. A complete analysis of the dynamics would 
require the numerical solution of the initial Gross-Pitaevskii equation. Such a numerical investigation of 
the problem should give definite answers to the following questions: 

(i) What are the conditions when the effect of the resonant excitation of coherent topological modes 
becomes possible? 

(ii) What types of coherent modes could be actually excited by this method? Could we create vortex 
modes in this way? 

(iii) What is the accuracy of the two-level approximation in the case of a resonant modulating field? 

(iv) What is an exact equation for the critical line, where dynamic critical phenomena may appear? 

(v) More generally, what are the critical manifolds in the parametric space, when there are several 
system parameters, and what is the relation of these manifolds to the bifurcation manifolds? 

Answering these questions would provide practical suggestions for creating various coherent modes of 
atom lasers. The corresponding numerical investigation is now in progress. 
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